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0> ■ Considered is quantum tunnelling in anisotropic spin systems in a mag- 

netic field perpendicular to the anisotropy axis. In the domain of small field 



< 

in 



the problem of calculating tunnelling splitting of energy levels is reduced to 
constructing the perturbation series with degeneracy, the order of degeneracy 
being proportional to a spin value. Partial summation of this series taking 
^3 . into account "dangereous terms" with small denominators is performed and 

o, 

the value of tunnelling splitting is calcaulated with allowance for the first 
. correction with respect to a magnetic field. 
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Tunnelling splitting of energy levels and the appearance of gaps m the energy spectrum 
of a quantum system is connected with breaking the symmetry of an unperturbed system by 



perturbation removing degeneracy. Usually, this effect is realized in quantum mechanics in 



two typical situations. The first one consists in that the Hamiltonian can be represented in 
the form H = H + V where H is the unperturbed Hamiltonian and V is perturbation, so 
correction to energy levels can be calculated with the help of some version of the perturbation 
theory. Another situation is tunnelling splitting of energy levels for a particle moving in a 
symmetric double well due to tunnelling penetration into a classically forbidden region. In 
that case the tunnelling splitting can be found by means of one of quasiclassical methods 
(WKB, instanton calculus, etc.). Although the physical reason - removing degeneracy - is 
the same in both cases, it is described in terms of essentially different languages regions of 
applicability of which, generally speaking, do not match. 
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It turns out, however, that there exists a whole class of physically interesting systems 
for which both approaches can be combined - anisotropic spin systems whose Hamiltonian 
is a quadratic-linear combination of spin operators. Their energy spectrum admits rigorous 
description in terms of an effective potential 0: energy levels of a spin system coincide with 
2S + 1 low-lying levels for a particle moving in a potential field of a certain form. If such 
a spin system is of the easy-type, for a sufficiently small field the corresponding potential 
has the form of a double well typical of tunnelling effects, so the energy splitting can be 
obtained, for example, with the help of instanton formulas fl[], H accuracy of which grows 
as a magnetic field decreases and the spin value increases. On the other hand, in the domain 
of small fields the perturbation theory can be applied. Thus, for sufficiently small magnetic 
fields and S ^> 1 both approaches match. 

The nontrivial peculiarly of the perturbation theory in such a situation consists in that 
the order of degeneracy is proportional to the spin value ||. As the quantity 5* plays the role 
of the inverse Planck constant h, it turns out that the degree of degeneracy tends to infinity 
in the quasiclassical limit. In so doing, the terms like B appear in the perturbation series 
which are nonanalytic with respect to S^in the limit S — > oo that in itself testifies that, 
although calculations are carried out on the basis of perturbation theory, the effect in fact 
admits the tunnel interpretation. 

In the leading order of perturbation theory the result of the calculation of tunnelling split- 
ting has been reported in ||] (without detailed derivation) and, independently, in where a 
consistent derivation of the formulas for energy splitting was suggested. This derivation was 
not direct application of standard formulas of the perturbation theory form textbooks and 
required careful singling out and summing all "dangerous" terms. In the present paper we 
perform the next step and calculate the first non- vanishing correction to the tunnelling split- 
ting for an arbitrary value of a spin. The approach developed in |I[ and the present paper is 
of systematic character that enables one to construct the series with respect to a magnetic 
field value including high-order corrections. In this sense the calculation of tunnelling rate 
is reduced to summing up the perturbation theory series with a given accuracy. 
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In recent years the interest to the nature of quantum tunnelling in spin systems of 
different nature sharply increased. It concerns the crossover between quantum and thermal 
domains |J, tunnelling in small ferromagnet particles |J, many-spin molecules |7[], etc. 
Therefore, formulae for the tunnelling splitting obtained below can be of interest not only 
from a viewpoint of general methods of quantum theory but for description of concrete 
physical systems observable in experiments. 

Let us consider the system described by the Hamiltonian 

H = -Si - BS X (1) 

Here Si are the operator of a corresponding spin component, B is a magnetic field. Since 
the first terms enters ([D with a negative sign, the system has an easy axis z. In the absence 
of a magnetic field all energy levels, except the one with S z = in the case of an integer 
S, are two-fold degenerate. A magnetic fields removes this degeneracy. The problem is to 
calculate the energy splitting due to quantum tunnelling with the first correction in B taken 
into account when B — > 0. 

Consider the Brillouin-Wigner perturbation series 

V 2 V 2 

E = e a + ^^ + ^^... (2) 

where e a = —a 2 corresponds to the unperturbed level with a ^-projection of spin equal to 
a, a a = E — s a , the perturbation V = —BS X and we took into account that diagonal matrix 
elements of perturbation are equal to zero. In order to determine AE a correctly, one needs 
to take into account explicitly the role of small denominators in "dangerous" terms of the 
perturbation series by performing a partial summation of (0). Following we sketch this 
procedure below. Consider, for simplicity, the splitting for the ground energy level AEq. 
The leading "dangerous" terms / reads 

r _ (Vs,S-lVs-l,S-2---V-S+l,-s) 2 /p s_i m\ 

T -[(E-es-i)(E-es- 2 )...{E-e-. s+1 )]* K 1 ' 

If we introduce a sequence of points on a number axis from S to — S corresponding to 
different values of a, then this term corresponds to a single sequential transition from S to 



— S and back. Now take into account that some elements of the number axis cab ne passed 
back and forward a few times (after a transition from S to —S the system returns to the 
intermediate point a, then to — S again and so on, passing the intermediate values of a). 
This suggests that additional factors appear in the perturbation series. Thus in the sum 
over terms which contain / as a common factor one should include higher orders of the 
perturbation. It is easy to guess that this leads to the following form of the last retained 
term in the perturbation series: 

t t ff-l i | 2 i \ f(-\ \-l (jp \-l V-s,s+\V-s+x,-s , 
f^f = f(t + r + r +...) = f[l-r) , r = (E - e s ) X, X = ~ + - 

ft — £-s+i 

(4) 

The factor x has the same form as the initial perturbation series @ that allows one to 
construct a simple equation with respect to E: 

(TP Vs,s-i Vs,s+i n2 2 Vs,s-iVs-i,s-2—V-s+i,-s 
(E - e a - — — — — h ...J = g s , g s ~ 



E — e S -i E — e s +i '' " (E — e S -i)(E — e s _ 2 )...(.E — £-5+1) 

(5) 

Generalization to other levels is straightforward and the result reads 

T/2 T/2 tt2<t-1t/ 

{E-e a -- + ...) -g a) g a - 2a-i n W 

ft — ^cr-l ft — £<r+\ ii fc=l a cr-fe 

The differences of the signs corresponds to the level splitting (the remaining terms coincide, 
making no contribution to the splitting). Taking advantage of the explicit form of the matrix 
elements V a)(T > we get 

AW ^=4> v ( ^'-i)- (7) 

in accordance with Q], ]I[ where a = S — n, n is the number of the unperturbed level. 

The first correction to (|7|) stems from three origins. First, the value of g a in (^|) acquires 
the correction due to field dependence of levels E a . Second, the correction appears due to 
E a which enter denominators in terms containing V 2 in the left hand side of (|6|). Third, the 
correction arises because of the next order terms in the perturbation series itself which change 



the structure of g a . Let us write down the total fractional correction as £ = = £i+£2+^3 
and calculate each term separately. 

The first term £1 = 1 ^l"^"^ ■ Substituting into this formula the explicit expres- 
sion for unperturbed levels = —a 2 and take into account the second order correction to 
E which can be easily found from the standard perturbation theory we obtain (n = S — a) 

B*[n*-2Sn + S(2S+l)] ^ 2S _ 2n _ lu _ 1 
41 " (25 - 2n)[(25 - 2n) 2 - lp =1 W 

The second type contribution to the correction can be found directly from the Brillouin- 

Wigner perturbation series by elementary methods: 

B 2 (n + l)(2S-n) n(2S - n + 1) 
^ 2 4 [ (2S-2n-l) 2 (25-2n + l) 2j l J 

The calculation of £3 is the nontrivial and rather cumbersome part of computations. This 

quantity is connected with the correction to the perturbation theory itself which is to be 

now cut off at terms of the order B 2S+2 . According to [0, a typical term of the perturbation 

theory series arises as follows. It is necessary to take the segment [— o~, a] and allow for 

all possible paths of the transition between these two points provided each an intermediate 

point on each path is being passed two times (to and fro). The main corrections to these 

terms arise if the transition from a given point to a next one contains one extra jump to 

and one - fro (delay and "marking time"). Summing up over all possible points - origins of 

extra jumps to and fro - we obtain 

& = efc -- fc+1 (10) 

Ota-kO-o-k+l 

Summing up all three contributions we find, after direct but somewhat lengthy calculations, 
the final expression which proves to be surprisingly simple: 

e = 1 - £ 2 7 , (11) 

= (2S + l) 2 (q+l) 
7 2(2a-l) 2 (2a + l) 2 

Thus, for an arbitrary energy gap (except the highest one in the case of semi-integer S) 

we get 



AF = (2g-n)! 2S , 2 _ (2S + inS-n + l) 

n 2 2S ~ 2n - 1 n\(2S - 2n - l)! 2 [ 2(25 - 2n - 1) 2 (2S - 2n + l) 2 J 1 ' 

This formula is supplemented for semi-integer S by the expression for the highest gap: 

AE s _ 1/2 = (S+ \)B[l - i(5 + - ^)B 2 ] (13) 

The obtained formulae, as was mentioned above, are valid when B — > 0. However, the 
presence of the correction term in (|T2"D and (|T3|) aids to make the area of applicability of 
these formulas more precise. First of all, the correction gives a qualitative notion about the 
approximation: it is seen that without account for such a correction the result is overstated. 
Further, the quantitative role of the correction turns out two-fold. 

First, it enables us to estimate the accuracy of the first approximation. Let us restrict 
ourselves by remarks concerning the ground gap (S > 1/2) 

AB » = ^sy B ^ - S^W B2] (14) 

It is seen from this formula that in the case S ^> 1 ,with the correction neglected, the 
fractional error is B 2 /8S. That means that even for B = y/S the computation of the ground 
gap gives the error approximately 10% while in the case B = 1 it proves to be very small, 
being of the order 0.1/ S. 

Second, the allowance for the correction essential increases the accuracy of computations. 
Thus, even for so large values of a magnetic field as y/S can be thought of to be, taking into 
account the correction enables one to make the fractional error one order smaller. 

It is worth noting that while computing energy gaps numerically even for moderately 
large values of S one is led to take into account very big number of digits that can be 
achieved only by modern computer computation systems. For example, in the case B — 1, 
S = 10 the gap has the order 10 -20 while for S — 50 it is of the order 10~ 200 . 

Either the main term or the correction in ( |14|) are in agreement with the instanton 
approach fl|] but for S ^> 1. Meanwhile, it is worth stressing that formulae (|T2"1) -(|13"1) are 
applicable for all values of S and embrace all the gaps. 
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In the case of an easy-plane spin system, i.e. a system with a Hamiltonian H = 
+Sf — BS X , the energy spectrum of which differs from that of the easy-axis type by the 
interchange of the sign of energy. Therefore, formulae for gaps remains valid provided they 
are renumbered, so that the lowest gap becomes the highest and vice versa. As far as the 
effective potential is concerned, in this case we have the model of a periodic potential energy 
in which the energy levels of a spin system coincide with the corresponding edges of energy 
bands (1). 

The elaborated approach describes quantum tunnelling in terms of the perturbation 
theory and enables one, in principle, to construct the energy level splitting as a series with 
respect to B with any given accuracy. 
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